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Abstract 

We study the spectrum of the QCD Dirac operator near zero virtuahty for Nc = 2. 



Qh! According to a universality argument, it can be described by a random matrix theory 

i-^ ' with the chiral structure of QCD, but with real matrix elements. 



Using results derived by Mehta and Mahoux and Nagao and Wadati, we are able 
obtain an analytical result for the microscopic spectral densi 
generating function for Leutwyler-Smilga type spectral sum rules. 



H I to obtain an analytical result for the microscopic spectral density that in turn is the 
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1. Introduction 

Quantum Chromodynamics (QCD) for two colors {Nc = 2) shares many common 
features with QCD for Nc = 3, generally believed to be the correct theory of strong 
interactions. However, there are important differences. In particular, the low-energy 
excitations do not only involve mesons but also baryons, which consist out of two quarks 
and are bosons. The corresponding effective theory [^ has a much richer structure with 
many Goldstone particles than for A^^ = 3. For example, for Nj = 2 we have nine 
Goldstone bosons |l| for the SU{A) — > 5*0(4) symmetry breaking scheme, instead of the 
usual 3 (see 0, ^ ^ for a detailed discussion of the different chiral symmetry breaking 
schemes for Nc = 2). 

As we have learnt from the work of Leutwyler and Smilga [^ , the existence of a low- 
energy effective theory imposes severe restrictions on the spectrum of the Dirac operator 
in the form of sum rules for the inverse powers of its eigenvalues. Recently, we have 
shown p, that the spectral sum rules can also be obtained from a random matrix 
theory with the chiral symmetry of the Dirac operator. This theory enabled us to derive 
the joint spectral density that generates all sum rules obtained by Leutwyler and Smilga. 
In fact, arguments from random matrix theory imply that the spectral correlations near 
zero virtuality are uniquely determined by the symmetries of the system 0, 0] . In other 
words they are universal (see |^ for a systematic study of this issue). This point has been 
known for quite some time in the study of spectra of classically chaotic systems |^, [10| [|11[] 
and the theory of S'-matrix fluctuations |T^. For example, both Ericson fluctuations [|I^ 
and universal conductance fluctuations [|14[ |15| can be unified within the latter context 
(see |jl6| for a discussion of this similarity). In particular, this means that the microscopic 
spectral density, defined as the V4 — i> oo limit of the spectral density while the eigenvalues 
are rescaled ~ V4, is a universal function. 

The random matrix theory that corresponds to the standard scheme of chiral symmetry 
breaking, has complex matrix elements. For this reason we have called it the chiral unitary 



ensemble (chGUE). In the framework of this model, new sum rules ||T^ can be derived 
with the help of the Selberg |T8|, |19[ integral formula. 



Because the effective theory for A^c = 2 is different from the generic case, we expect 



a different microscopic spectral density and different sum rules. This raises the question 
what is the correct random matrix theory for A^^ = 2. The answer becomes clear if one 
considers the matrix elements of the Dirac operator. In this case it is possible to choose 
a basis in which they are real, as opposed to three and more colors where they are com- 
plex. This reminds us of the three universality classes in random matrix theory pO[, the 
Gaussian Orthogonal Ensemble (GOE), the Gaussian Unitary Ensemble (GUE) and the 
Gaussian Symplectic Ensemble (GSE). The first correspond the real symmetric matrices, 
the second to Hermitean complex matrices and the third to quaternion matrices. It is clear 
that the correct random matrix ensemble should not only embody the chiral symmetry 
of the Dirac operator but also satisfy the additional constraint that the matrix elements 
are real. From now on we will call this matrix ensemble the chiral orthogonal ensemble 



(chGOE). The third possibility is realized for fermions in the adjoint representation |21 
In this case it is possible to regroup the matrix elements of the Dirac operator in terms 
of quaternions. The corresponding ensemble will called the chiral symplectic ensemble 
(chGSE). 



In a separate publication |2^ we will show that the cliGOE describes the spectrum of 



the Dirac operator in a liquid of instantons [^, whereas it is not given by the chGUE 



In this model the Dirac operator is diagonalized in the space of fermionic zero modes. 
Indeed, the matrix elements of the Dirac operator are real. Also the corresponding sum 



rules will be given elsewhere [17 



In this work we report on the calculation of the microscopic spectral density (defined 
in section 2) of the Dirac operator for Nc = 2. The argument that the matrix elements of 
the Dirac operator are real is presented in section 3. The corresponding random matrix 
theory and the joint eigenvalue distribution is derived in section 4. The exact spectral 
density for any finite size random matrix is calculated in section 5. The microscopic 
spectral density is obtained in section 6 and concluding remarks are made in section 7. 
Some technical details are worked out in three appendices. 



2. Formulation of the problem 



The distribution of the eigenvalues of the Dirac operator are determined by the fluc- 
tuations of the gauge field which are subject to the Euclidean QCD partition function 

Nf 

z=j:e-'{ii n(A^+^?K)^.(^)' (2-1) 

•^ /=1A„>0 

where the average (■ ■ ■)s^{A) is over gauge field configurations with topological quantum 

number i/ weighted by the gauge field action Si,{A). The product is over all eigenvalues of 

the Dirac operator, and relevant observables are obtained by differentiation with respect 

to the masses rrif. The number of flavors is denoted by Nj. The factor expiOu represents 

the topological term in the action. 

The condensate can be expressed as a derivative of the partition function 

i d 
(gg) = lim lim — \ogZ{mf). (2.2) 



According to the Banks-Casher [^ formula we have 



m = -^. (2.3) 

where the spectral density p(A) is defined as 

p{X)=J26{X-Xn). (2.4) 

It is now clear that, in order to obtain a nonzero value of {qq)j we should have 

(p(0)) ~ \/4, (2.5) 

or, put differently, the spacing between the eigenvalues near zero virtuality is ~ I/V4 
(Note that for a non-interacting system the spacing between the eigenvalues is ~ 1/V^^^). 
As was observed by Leutwyler and Smilga |0|, (2.5) implies the existence of a family 
of new sum rules. The simplest one involves the sum 



^4 A„>0 \'^n/ ^ 



which should converge to a finite limit for V4 — » cxd. 

The above mentioned sum-rules can be expressed in the microscopic spectral density 
defined by 

Psix) = hm i (p(^ )\ (2.7) 



in the sector of topological charge v. For the sum (2.6) we find 

/dx^. (2.8) 

In this paper we will obtain an analytical expression for ps{x) for Nf flavors and arbitrary 
topological charge u. 



3. Symmetries of the Dirac operator for Nc = 2 

In this paper we study the Euclidean Dirac operator 

D = ijd + jA, (3.1) 

where A is an SU{2) valued gauge field. The spectrum of D is defined by the eigenvalue 
equation 

D(Px = A0A. (3.2) 

The Dirac operator for SU{2) has two symmetries. First, the chiral symmetry, which is 
present for any S'f/(A^c)-valued gauge group, 

{l5,D} = 0, (3.3) 

and second, a symmetry specific to SU(2): 

[C-'r2K,D] = 0, (3.4) 

where C is the charge conjugation operator (7* = —C^muC^^ and CC* = —1), and K 
denotes the complex conjugation operator. Because of the first symmetry, the eigenvalues 
occur in pairs ±A. The second symmetry operator has the property that 

iC-'T2Kf = 1. (3.5) 

As is well known from the analysis of the of the time-reversal operator in random matrix 
theory (see for example refs. pSJ |T9l), this property allows us to choose a basis in which 



the Dirac operator is real. In a chiral basis the Dirac operator therefore has the following 
general structure 

(11)^ (3.a) 

where T is a real matrix. This is the only information of the Dirac operator that will be 
injected into the random matrix model to be defined in section 4. 



4. The chiral random matrix model 

As has been shown in previous section, the matrix elements of the Dirac operator can 
be chosen real. The corresponding random matrix theory with the chiral structure of 
QCD is defined by the partition function 

Nf 



Z„^/l>TP(r)lidet(;:V ^;^J, ,4.1) 



/ 

This model is defined for Nf flavors with masses rrif in the chiral limit {nif —>■ 0) and in 
the sector with topological charge z/. The latter property is implemented by choosing T 
anm X n matrix (for definiteness, m > n) with \m — n\ = u. With this choice the matrix 
of which the determinant is calculated has exactly z/ zero eigenvalues for ruf = 0. The 
integral is over all matrix elements of T, i.e., VT is the Haar measure. As follows form 



the the maximum entropy principle |^ the distribution function of the overlap matrix 
elements P{T) is chosen Gaussian 



2^ 



P(T) = exp(-^TrTf). (4.2) 



In the present case (T real) the value oi (3 = 1, whereas for complex matrix elements 
(3 = 2. With this convention, a mean-field argument shows that the average level density 
for the chGOE is the same as for the chGUE (see appendix A). One can easily show that 
the eigenvalues of a matrix with this block structure occur in pairs ±A, a property that 
is well-known for the nonzero eigenvalues of the massless Dirac operator. The density of 
modes N/V^ {N = 171 + 11) is taken equal to one, which allows us to identify A^ with the 
volume of space time V4, which we will do from now on. In agreement with general QCD 
relations we will always assume that u <^ N. 



The matrix T can be diagonalized by an n x n orthogonal matrix Oi and an m x m 
orthogonal matrix 02'- 

T = O1AO2. (4.3) 

Here, A is an nxm diagonal matrix with diagonal matrix elements A^. The joint eigenvalue 
distribution is obtained immediately by using A, Oi and O2 as new integration variables. 
Because the integrand only depends on A, the integration over the orthogonal matrices 
can be performed trivially. The Jacobian corresponding to the transformation (4.3) is 
given by (see appendix B) 

n n 

J{^) = \1\>H.->h\\IK (4.4) 

k<l k=l 

The joint probability density of the nonzero eigenvalues is therefore given by 



n 



p(Ai,---,A„) = niA^-AfinAr^''^exp(-^EAD- (4.5) 

k<l k 

Since this distribution is symmetric in all eigenvalues, the spectral density is simply ob- 
tained by integrating it over all eigenvalues except for one: 

J rfAi---rfA„p(Ai,---,A„) 
The normalization integral in the denominator will be denoted by Z{n). 



5. Calculation of the spectral density 

To make contact with what is also called the orthogonal generalized Laguerre ensemble, 
we introduce new integration variables by 

Xk = — |. (5.1) 

If we absorb the constants in a redefinition of the normalization constant Z, the spectral 
density is given by 

p{\i)d\i = pL{xi)dxi = —- j dx2- ■ ■dxnW\xk-xi\Wxle:!Ci>{-^Xk)- (5.2) 

^ •'^ k<l fc=l fc=l 



For convenience we have introduce the notation 

a = Nf-\ + '^. (5.3) 

General expressions for spectral densities given by these types of integrals have been 
derived by Mahoux and Mehta p^. The result is expressed in skew-orthogonal polyno- 



mials Rk defined by 

{R2k,R2i) = {R2k+i, R21+1) = 0, (5.4) 

with the scalar product given by 

/•OO /"OO 

{f,g)= dxx^'e-^ dyy^e-ye{x-y)f{x)giy), (5.5) 

Jo Jo 

where e{x — y) = |sign(x — y). For this weight function the normahzation constants r^ 

are known explicitly 






-^ = htnt' = ^ ^5^W^ '-, (5.6) 

where hl"^^ are the normalization constants of the monic generalized Laguerre polyno- 
mials with index 2a + 1. They are fixed by the normalization integrals Z{n) (see below 
(4.6)). 

The result for the spectral density is given by 

(n/2)-l ^ 

pLix)= Y. — (02m(a;)02m+i(a;)-02m+i(a;)02m(a;)), (5-7) 

where 



m=0 ^rn 



("CXD 

M^)= dyy''e-ye{x-y)Ru{y). (5.8) 

Jo 

From the definition of the skew-symmetric scalar product one finds for the normalization 

I^ PL{x)dx = n. 

The polynomials Rk can be obtained most conveniently by expanding them in the 

monic generalized Laguerre polynomials C'^{x) = n\L'^{2x)/{—2)": 

2m 

R2m{x) = 2^ amfcC'2m-fc(^)) 
fc=0 
2m+l 

R2m+lix) = J2 KkCt,^,.k{^)- (5.9) 

fc=0 



Starting from the identity 



oo 



((- - 1)/, g) - (/', g) = - x''^e-'^fix)gix)dx (5.10) 

X Jo 



Nagao and Wadati I^Q] where able to derive recursion relations for the expansion coeffi- 
cients 

Clmn+l = --{'2m-n)amn, U > 0, 

bmn+1 = --{2m-n + l)bmn, n>2, 

in 
bm2 = - — {2bmi + 2a + 2m + l), (5.11) 

where bmi is not determined by the recursion relation, and in fact, does not contribution 
the spectral correlation functions. The intial conditions are fixed by 

O-mO = bmO = 1- (5.12) 

The solution of these recursion relations is straightforward: 

1 f 2m' 



(-2)" \ n 

2a + 2m + l f 2m 



n\, (5.13) 

n-1)!, n>2 (5.14) 



(-2)"-! \n-l ^ 
where we made the choice (does not satisfy (5.14) for n = 1) 

6^i = -(2a + 2m + l). (5.15) 

Remarkably, using the identity Z^fcLo-^m-fc = -^m"^ (^^e [Q) the sums in eq. (5.9) can 
be performed exactly. The result is 



2myX) — i^2m 

2a + 2m + l 
2 



R2m+i{x) = Ct^,{x)+ ^, ^ iCl:^^\x) - mCtt\{x)) (5.16) 



The spectral density (5.7) can be written as the sum of two terms 

^{n/2)-l 



p,{x) = x'^e-^ dyy'^e-ye{x~y)l ^ -iCt\\i^)Crn^\y) ^ Ct^\^)CtX\iy)) 

m=0 ^'"™ / 

(5.17) 



The two sums can be combined into a single sum. This expression can be simphfied 
further by applying the inverse Christoffel-Darboux formula (see ||3^) and collecting all 
terms. We find 

Pl{x) = x^e-^ ry^e-y\^ £ ^^^^g^Cr^(x)Cr^(y) (5.18) 

In order to take the limit n ^ oo the n-dependence has to be made more explicit. To 
achieve this, we apply the Christoffel-Darboux formula once again to the term proportional 
to n, and, after using the identity, 

m + lax 
to the term proportional to [m + 1). The term proportional to n cancels against one of 
the terms obtained from the differentiation with respect to x and y. We finally obtain 

Pl{x) = x^e / dyy^e ^6 x - y) [ ^f^ ^ ^^^L^^ 

r(n + 2a) Jo \ [x — y)^ 

^ Ll^-\2x)Ll^+\2y) + Ll^^\2x)Ll^-\2y) - 2Ll%2x)Ll^{2y) \ _ 

x-y )' 

(5.20) 



This sum can be written more compactly as an integral over the unitary kernel K{2x, 2y) 

d d 
'dy 

with the kernel K{2x, 2y) defined by 



f°° a a 

Pl{x) = x^e-- j^ dyy^e-ye{x-y){--—)K{2x,2y), (5.21) 



TTfo 2^-n! Ll%{2x)Ll^{2y)-Ll%2x)Ll%{2y) 

A(2x,2yj-^^^^2«) 2x-2y ' ^^'^^^ 



This kernel was first considered by Fox and Kahn |3^, and was studied in great detail 



by Bronk ||3T|. The relation (5.21) shows an intimate and, an as yet not well understood. 



relation between the unitary kernel and the orthogonal spectral density. 



6. The microscopic limit 

In this section we derive the microscopic limit of the spectral density. First, we express 
the parameter a in the mean level density. An expression suitable for the analysis of the 

10 



spectral density many level spacings away from the origin but yet far from the edge of 
the semi-circle is obtained by commuting the e-function through the derivative operators. 
The result can then be written as the sum of the chGUE spectral density plus a remaining 
oscillatory term: 

Pl{x) = 2x^''e-^''K{2x, 2x) + x^e"" / dyy^e-'^ij- - —)e{x - y)K{2x, 2y). (6.1) 

Because 

roo 

2 / dxx'^''e-'^''K{2x, 2x) = n, (6.2) 

Jo 

the second term does not contribute to the total number of levels. As is well known for 
the chGUE the large n-limit of the first term is a semicircle. In the normalization (4.2) of 
the distribution of the matrix elements, the average level density does not depend on j3 
(see appendix A). The second term in (6.1) therefore does not contribute to the average 
level density. From the asymptotic formula for the Laguerre polynomials (see eq. (C.2)), 
it follows that in the thermodynamic limit {n ^ 1) 

Pl{x) ^==, (6.3) 

vr \J 2nx 

for a; — !► 0. The large-n limit of the spectral density (4.6) at A = is then given by 

p(A = 0) = pz.(x)^ = — . (6.4) 



According to Banks-Casher formula |2J], we can identify the parameter a as 



a=l (6.5) 

where S is the chiral condensate. 

The microscopic limit of the spectral density is given by 

z d\ z dx 

ps{z) = lim p{X = —) — = lim pl{x = - — -) —. (6.6) 

N-*co jy dz N^oo Sna^ dz 

In order to evaluate this limit we start from expression (5.20) for the spectral density and 
use the microscopic variables z and w defined by 

z"^ w'^ 

X = ^ — ^, and y = - — -, (6.7) 

11 



where z is related to the original eigenvalues by z = AA^. 

The term in the integral proportional to L^"+^L^"~^ does not satisfy the conditions 
necessary to interchange the limit and the integration. However, if we add and subtract 
the term 



roo 



T{n + 2ay 

to the integral, it can be proved by dominated convergence that in the subtracted integral 
the limit n ^ oo can be taken before integration. The integral (6.8) has to be performed 
exactly first. The result is (see appendix C) 

/•OO 

Jimn-'^+W dyy''-^e-yLl^+\2y)=2~''+\ (6.9) 

From the asymptotic formula for the Laguerre polynomials (C.3) it then follows that the 
microscopic limit of the spectral density is given by 

wJ2aizJ:)J2a-liwJ:) - Z Jsa-l (^S) J2a(wS) 
^ (Z^)2«(z2 _ ^2) • y'^- UJ 

Note that half of the subtracted term is again reabsorbed into the integral. To achieve 
this we have used the identity 



-J2.+l(zS) = -- / dwJ2a^^{wT.)z'^ — {z~^\J2a{^z)). (6.11) 

4 Jo dz 



In eq (6.10) we have expressed the microscopic spectral density in terms of an integral 
over the Bessel kernel. This kernel was studied extensively by Widom and Tracy [p^ . 

The leading term in the small z-expansion is obtained by replacing z^ — w"^ ^ —w"^ 
and e{z — w) — > — |. Then all integrals can be performed exactly, and after a cancellation 
only the integral (6.11) and the first term in (6.10) contribute to leading order. This 
results in 

p(^) ^ -(—)"^''\ (6.12) 

^^ ' 2r(2a + 2) \2 J ^ ' 

The simplest spectral sum rule is given by 

r ^ Ps{z) S^ 

L ^^—-siNiT^ry ("3) 

12 



This sum rule can also be derived from the partition function without reference to the 
spectral density 



In the case of zero flavors {Nf = 0) and zero topological charge (z/ = 0) a direct nu- 
merical simulation of the random matrix model (4.1) is particularly simple. To convince 
the reader that (6.10) is correct we show in Fig. 1 a histogram obtained from the diag- 
onalization of 10,000 128 x 128 matrices (dashed curve) and the exact microscopic limit 
(6.10) for a = (full line). Perfect agreement is no exaggeration in this case. In Fig. 2 
we show the results of the microscopic spectral density for Nf = 1,2 and 3. We observe 
much less oscillations than in the chGUE case. This agrees with general property known 
from the classical random matrix ensembles that spectra of real matrices are much less 
rigid than spectra of complex matrices. Therefore, the variation of each level about its 
average position, called level motion, is much larger resulting in the (almost) absence of 
oscillations. 



7. Conclusions 

In this work we have studied the spectrum of the QCD Dirac operator for an SU{2)- 
valued gauge field in the fundamental representation. This case differs from any other 
nonabelian gauge group in the fundamental representation by the reality constraints of 
the eigenfunctions: it is possible to choose a basis in which the matrix elements of the 
Dirac operator are real. As we have argued before for the chGUE case, the correlations 
between the eigenvalues of the QCD Dirac operator near zero virtuality are insensitive 
to the detailed dynamics of the system and can be described by a random matrix model 
with as only its symmetries as input. There is no reason to believe that the present case 
is different. However, the appropriate random matrix ensemble not only has the chiral 
structure of QCD but also has real matrix elements. For this reason, it will be called the 
chiral orthogonal ensemble, abbreviated by chGOE. 

Because of the spontaneous breaking of chiral symmetry, the spectral density near 
zero is ~ V4. This property allows us to define a limit V4 — » 00 of the spectral density in 
which the eigenvalues are at the same time rescaled by a factor V4. The resulting spectral 

13 



density, called the microscopic spectral density ps{z), is insensitive to the dynamics of 
the system. It is a universal function that is entirely determined by the symmetries of 
the Dirac operator. 

In this paper we have evaluated ps{z) for the case of SU{2) in which the Dirac operator 
is real. This turned out to be much more difficult than for a complex Dirac operator, a 
well known property of random matrix theory. The spectral density in the present case 
differs from the chGUE case by the absence of strong oscillations. As is also the case for 
the classical random matrix ensembles, the spectrum of a complex matrix is much more 
rigid than the spectrum of a real matrix. This can be made more quantitative in terms 
of the so called level motion which turns out to be much larger for the chGOE than for 
the chGUE. 

The microscopic spectral density is generating function for the Leutwyler-Smilga sum 
rules. The sum rules for for Nc = 2 have not been obtained before. In view of the fact that 
the corresponding static effective field theory involves both baryons (which are bosons in 
this case) and mesons |^, it is not surprising that the results are different from for other 
nonabelian gauge groups in the fundamental representation. It would be interesting to 
derive the sum rules from the static limit of the effective field theory also in this case. In 
the case of one fiavor the effective theory for Nc = 2 and other nonabelian gauge groups 
coincides. Although the spectral density is different, all spectral sum rules for Nf = 1 
turn out to be the same for the chGOE and chGUE. Further work to clarify this issue is 
in progress 



Appendix A 

The joint eigenvalue density valid for both the chGOE {f3 = 1) and the chGUE (/? = 2) 
is given by 

p,(Ai,---,A„) = C,,.ni^'-^'l'nA^xp(-^EAD. (A.l) 

where Cn,i3 are normalization constants and a = {2Nf + (31/ + (3 — 1) . 

14 



The normalization integral / dXi ■ ■ ■ d\nPp{\i, ■ ■ ■ , A„) can be approximated by 

n/5 



exp 



(3 J dXdX' log \X^ - X'^\p{X)p{X') + J dXp{X) l-^X^ + alogX + pp{X) 



(A.2) 



where the average level density p{X) satisfies a 'mean-field' equation obtained by mini- 
mizing the exponent: 

2(3 [ dX' \og\X^ ~ X'^\p{X') - ^X^ + ahgX + p = 0. (A.3) 

The normalization of the level density is introduced via a Lagrange multiplier. By differ- 
entiation with respect to A^ we obtain the principal value equation 

r dX' , , „ nd a , , ■, 

Since J p{X)dX = n the third term is subleading for n ^ oo (in agreement with general 
properties of topological fluctuations we have v <^ n). Consequently, the 'mean-fleld' 
equation for p does not depend on 13. A more elaborate discussion of this argument for 



the classical random matrix ensembles can be found in the book by Mehta ||T^ 



Appendix B 

In this appendix we calculate the Jacobian of the transformation of the matrix valued 



variables T into its eigenvalues and eigenangles using techniques developed in |5^. For 
an arbitrary real n x m matrix we have 

T = O1AO2, (B.l) 

where A is a positive n x m diagonal matrix and, the n x n matrix Oi and the m x m 
matrix O2 are orthogonal (for deflniteness m > n). By differentiation of (B.l) it can 
be shown that that the variation dT can be expressed in variations 60i of Oi near the 
identity 

Cr^rfTOa = 60iA + dA-A602, (B.2) 

O^^dfOi = -A60i + dA + 602A, (B.3) 

15 



where we have introduced 50i = O^^dOi. Note that the matrices 50i are anti-syminetric, 
and, in particular the diagonal matrix elements are zero. For the invariant distance we 
find 

TrdTdf = Y^idAkf + E o [(^^1 - ^(^2)kif (A, + \if + Y.1; [(^^1 + ^02)Hf (Afc - kf 

k k<l ^ k<l ^ 

n m 

+ E E [(^0i)d'A^ (B.4) 

k=ll=n+l 

This allows us to read off the Jacobian of the transformation to the variables dKk and 
{50i ± bO'i)l \/2 from the Lame-coefficients: 

n n 

k<l k 

Note that the total powers of lambda can be verified by a dimensional argument. 



Appendix C 

In this appendix we consider the large n limit of integrals of the type 

/■oo 

ja,k^^k+i^a Ll^-'^^2y)y'^e-ydy, (C.l) 

Jo 

where k is an integer, and n is even. From the asymptotic expansion for n ^ oo of the 



generalized Laguerre polynomials (see ||3^ for this and other properties of the generalized 



Laguerre polynomials used below), among others. 



Vz(x)r^—^e'^x 2 47^2 4 cos(2v^x ), (C.2) 

n\ , r^ \ y 4 4 



it is clear that for k < —3/2 and a > — 1 the absolute value of the integrand can be majo- 
rated by an integrable function. In that case we can, after introducing a new integration 
variable hj y = w'^/2n, interchange the limit n —>■ oo and the integration over w. Using 
the asymptotic result 

lim L"(— ) ~ n''w-''JJ2w), (C.3) 

n^oo fl 

with Ja a Bessel function, we find that 

lim I^^" = 2-" / dww''+'J2a-k{2w). (C.4) 
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This integral can be evaluated analytically, resulting in 



,a,k o-a-ir(a + l) 



lim /„"''= = 2-°-^ ^ , ■ (C.5) 

r{a-k) ^ ' 



n— >oo 



For — I < A; < — I (since k is an integer: k = —1) and a > —1 the integral still 
converges but no longer satisfies the conditions necessary to interchange the limit and 
the integral. In this case there is an important contribution from the region around the 
largest zero of L". However, for /c = —1, it is particularly simple to evaluate the integral 
for any finite value of n. From the recursion relation 

n 

Ll-+\2y) = Y: L'^iy)L'^_„Xy) (C.6) 

we can reduce the integral to a normalization integral for the Laguerre polynomials re- 
sulting in 

L -n ^^^^^^ . (C.7) 

The asymptotic limit follows immediately from Stirlings formula 

lim P'-i = 2-\ (C.8) 

which is a factor 2 bigger than the result given in (C.5). We conclude that half of the 
contribution to this integral is from the region near the largest zero of the Laguerre 
polynomial. 
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Figure Captions 

Fig. 1. The microscopic spectral density ps{z) versus z for A'^^ = 0, z/ = and S = 2. The 
full line represents the exact analytical result and the cashed curve shows data obtained 
by diagonalizing 10,000 random matrices distributed according to (4.1). 

Fig. 2. The microscopic spectral density ps{z) versus z for Nf = 1 (dotted curve), Nf = 2 
(dashed curve) and Nf = 3 (full curve) all for z/ = and S = 2. 
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